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I •        IIIIOD UCTI C N 

Accurate  and  reliabla  rumerical  weather  predic-.ior.  is  a 
n.ajcr  goal  fcr  civilian  and  mili-ary  a- mospheric  studies. 
V9rifica-icn  studies  ar€  an  imfortan*  aspect,  cf  numerical 
model  d9v<=lcpirer:t  as  they  isolate  model  strengths  and  weak- 
nesses and  identify  errcrs  which  must  ba  eliminated  tc 
extend  the    present  limits    cf   forecast   skill. 

Ihis  thesis  will  investigate  the  Naval  Operational 
Global  Atmospheric  Prediction  System  (NOGAPS)  mcdsl's  abil- 
ity tc  forecast  the  pcsiticn  and  movement  of  SOOmb  waves  in 
".he  irid  latitudes  of  the  Northern  Hemisphere.  NOG  APS,  which 
becaie  operational  in  Sept=mb3r,  1982  at  the  Fleet  Numerical 
Oceancgraphy  Center  (FNOC)  has  b<=en  a  project  of  the  Naval 
Elnvircnmental  Prediction  Hesaarch  Facility  (NEPRF)  since 
•^976.  Th€  long  term  goal  of  NOGAPS  is  to  achieve  a  medium 
range  (7-1C  days)  numerical  forecast  capability  for  the  U.S. 
Navy  (F.osmcnd  ,  198 1)  .  Verification  cf  NOGAPS  is  very  impor- 
tant at  this  stace  due  tc  the  research  type  background  of 
thr  mcdel.  The  NOGAPS  fcrecast  model  is  a  modified  form  of 
the  CCLA  general  circulation  model.  Consequently,  it  has 
not  had  the  tenefit  of  irany  prior  intensive  verification 
studies  as  have  most  second  and  third  generation  operational 
models. 
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Eecsnt  studies  have  raised  serious  questions  conc^rr.ing 
the  relative  accuracy  of  the  planetary-scale  (zonal  wavenuin- 
bers  1-2)  viav€  mction  forecasts  at  500mb  produced  by  nunieri- 
cal  weather  prediction  models.  Lambert  and  Meril^es  (1978) 
corcluded  the  planetary-scale  motion  prediciions  of  a  spec- 
tral numerical  mcdel  were  poorer  with  respect  -o  persistBnce 
xhan  long  ( wavsrumbers  3-5)  or  medium- scale  (wavenumbers 
6-10)  motions.  Miyakoda  et  al.  (1972)  examined  the  GFDL 
model  perfcrmence  over  a  series  of  win-2r  forecasts  ani 
found  similar  error  characteristics  in  -he  planetary-scale 
forecasts.  This  study  will  verify  NOGAPS  planetary  wavs 
forecasts    and    other    wave    regimes,    large   and    medium-scale. 

500mb  data  from  toth  the  forecas*  and  the  analysis  are 
spirctrally  deccmpcsed  intc  i*3  basic  ccmponen-s  of  ampli-.  ude 
and  phase  then  organized  in-o  planetary,  long  and  medium- 
seal?  waves  for  separate  analysis  of  their  atmospheric 
be'iavioral  characteristics.  If  one  assumes  spatial  coupling 
(the  interacting  cf  different  scales) ,  then  mcdel  weaknesses 
in  forecasting  'long*  wave  motion  will  be  reflected  in  the 
accuracy  of  'short'  wave  motion  forecasts.  The  purpcss  cf 
this  study  is  tc  perform  a  thorough  analysis  cf  the  fcr^- 
CaS-t£  of    each    wave  grouping.      Somerville     (1980)     examined    th^ 
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plan€tary-scal€  wave  forecasts  cf  a  primitive  equation 
numerical  model  in  both  a  global  and  hemispheric  configura- 
tion. He  concluded  that  the  global  vsrsion  was  aark^dly 
more  skillful  than  th€  hemispheric  version  especially  in  ths 
lat-€r  part  cf  the  five-day  forecast  period.  While  the 
present  study  will  consider  only  a  glcoal  version  of  NOGAES 
tne  results  will  lend  themselves  to  an  application  of  Scoier- 
ville's    work    to    NCGAPS. 

E::rors  in  phase  speeds,  amplitudes,  development  of  bare- 
clinic  systems,  damping  of  smaller  scale  features  and  plane- 
tary wave  structure  become  clearer  when  the  forecasted 
variarles  are  analyzed  over  time  via  a  Hovmoller  or  trough- 
ridge  plot.  Baumhefner  and  Downey  (1978)  utilized  -his 
technique  tc  interccmpare  the  forecasting  skills  of  three 
numerical  weather  prediction  models.  They  constructed  lon- 
gitude-time plots  along  a  fixed  latitude  to  better  under- 
stand the  tine  evclution  cf  the  forecast  fields.  Forecasts 
waich  originally  appeared  quite  good  were  revealed  tc  have 
amplification  errors  in  many  of  the  transient  systems.  A 
spectral  tr=iakdcwn  cf  the  Hovmcllar  plot  in  planetary 
(wavenumbers  1-3),  large  (wavenumbers  5-7)  and  medium 
scale-waves     (wavenumbers    8-12)    were    utilized    in    their    study. 
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On*    has   only      to    consider    the   complexity    of      a   typical    SOCtnb 
wave-    pattern    tc    appreciate   zhe   benefits  of    this   approach. 
The   specific   cbjectives  cf  this   thesis   are: 

(1)  E'e  velopment  of  a  spectral  verification  program  using 
wavenuinber  grcuping  and  Hovmoller  plots  for  trcpospheric 
height    data, 

(2)  Application  of  the  program  to  interesting  NOGAPS  fore- 
cast   sit.uaticns    at   SOOmb, 

(3)  Use  of  wavenumber  analysis  to  verify  NOGAPS  forecasts 
of    plan-5tary,    long  and   medium-scale    waves. 

Chapter  2  reviews  the  main  features  of  the  NOGAPS  model. 
The  pafiern  decomposition  program  is  described  and  demon- 
strated in  chapter  3.  Chapter  4  presents  the  basic  charac- 
tsri sties  cf  the  case  study  periods  and  their  subsequent 
atmospheric    behavior.  finally,      chapter    5   •  summarizes    the 

results    of    the    forecast    ccirpa  risen   and    makes    suggestions    for 
further    research. 
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II.        NCGAPS    SUMMARY 

Iha  Naval  Cp<=ratianal  Global  Atmospheric  Predic-ion  Sys- 
tem (SOGAFo)  was  developed  by  the  Naval  Environman-al  Pre- 
dict icr^  Research  Facility  (NEPRF)  to  be  run  at  the  Flee- 
Numerical  Cc^ar.ccraphy  Cer^.er  (7N0C).  The  objectivsrS  of  the 
NOGArS  project  are  to  ianprove  tropical  and  Southern  Hemis- 
peric  forecasting  and  achieve  mid-range  (5  -  10  days)  global 
numerical  forecast  capability  to  support  U.S.  Naval  opera- 
tions    (Rosniond,1981)  . 

The  objective  analysis  method  used  by  NOGAFS  is  a  varia- 
tion cf  the  successive  corrections  technique  {Barnes, 1954) 
which  discr iiinates  among  observing  systems  according  to 
their  accuracy.  While  the  final  analyses  is  net  as  scchis- 
tic5t€d  as  the  octimal  interpolation  method  it  requires  far 
less   ccmputer    resources. 

GloDal  mass  and  wind  fields  are  analyzed  independently 
with  cpticnal  data  bogusing.  Eata  rejection  criteria  are 
intentionally  Icose  since  -he  initialization  scheme 
(Barker, 1 981)  removes  the  gravity  waves  excited  by  bad  data. 
This  technique  prevents  inadvertent  rejection  of  the  good 
data    associated    with    rapidly    developing   systems. 
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A  calculus  cf  variaticn  methcd  combinss  th^  rasults  of 
the    Viind    and    mass  analyses.  The   hydros'aiic   equation    per- 

niits  sabstituticn  cf  teufsrature  for  gsoporr-n-ial  which 
gives  the  forecast  model  three  dimensicnality .  The  form  of 
NOGAES  init:.alizatioD  used  in  our  study  produces  balancsd 
geopct£-u-ia:.s  (temperatures)  and  non-di  vsrg  =nt  winds  which 
are  in--3r pc.'.ated  directly  ro  the  model's  sigma  coordir.a.e 
system.  i:ie  method  which  became  operational  in  December, 
1982,  defi:ies  geopotenrial  and  wind  as  objective  analyses 
corr-c'.ion  fields.  The  resul-ing  balanced  corrections  ars 
in-^-erpolat  e  3  tc  -he  sigma  surfaces,  as  in  rhe  firsx.  method, 
ar  d  ?dc'}d  tc  the  first  guess  forecast..  This  method  pr  =  - 
servei?  th9  model-generated  divergence  and  minimizes  the  ver- 
tical  int  er  pclat  icn    error. 

Ihe  first  guess  forecast  of  specific  humidity  and  the 
prcancitic  plenetary  boundary  layer  (?BL)  variables  are  used 
as  iritial  values  for  the  subsequent  forecast,  cycle.  The 
six-hcur  update  cycle  is  preferred  over  -he  twelve-hour  con- 
figurstion  since  Isss  tine  inter pcla tion  of  asynoptic  data 
is  required.  However,  the  shorter  cycle  is  more  subject  to 
initialization  shock  and  relies  heavily  on  the  initializa- 
tion   technique    to   control    the    noise.  Failure    to    do    so    may 
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rssui"^.  in  suts€qu€nt  first  guess  errors  suffici-ntly  lirge 
to   cause    rejecticn  of   good  observational   da-a. 

Ihe  NOGAES  model  is  a  modified  version  of  the  UCLA  gen- 
eral circulcticn  model  (GCM)  described  oy  Arakawa  and  Lamb 
(1977).  Ihe  2.4  by  3.0  degr9=?  horizon-al  grid  uses  scheme  C 
s-agg€ri:.g  of  the  variables  and  conserves  energy  and  enstro- 
ph y  khen  the  flow  is  condivergen t .  All  dependent  variables 
except  v-~-rtical  velocity  are  defined  in  the  aiddle  of  each 
of    the   six    vertical    sigma    levels     (Fig.    1)  . 

Moisture,  racist  static  energy  and  moaientun  are 
assumed  to  be  well-mixed  in  th=:  PEL.  The  aquations  govern- 
ing the  diatctic  processes  are  based  en  the  method  cf  Dear- 
dcrff  (1972).  These  equations  represent  the  conservation  cf 
mass  in  -.h€  ?EL,  temperature  and  specific  humidity  inversion 
streigtb  ar.d  the  jump  in  momentum  across  the  inversion.  A 
special  feature  of  the  model  is  the  variaole  depth  of  the 
?BL  in  the  lowest  layer  which  has  a  capping  inversion  cf  a 
porous      .material    surface.  This   allows      entrained   mass      -^o 

deepen  the  PcL  and  removed  mass  (by  cumulus  mass  flux)  to 
decrease  the  FBI.  The  fluxes  at  the  top  of  the  PEL  are  net 
a  constant  fraction  of  the  surface  fluxes  but  rathar  a  func- 
tion     of    the      entrainment .  Dissipatina      stratus   causes      a 
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Figura    1.       The    Vertical    Differencing   of  the    NOGAPS    Model 


17 


par-.iculaily  large  flux  across  -he  inversion  fcr  sira^us 
capped   PBl's. 

a  significant  riodif ication  re  the  GCM  is  -he  cons-rain- 
ment  cf  -he  FEI  -c  remain  in  -he  bottom  sigma  lay=r.  This 
eliminates  serious  computational  problems  yer  has  no-  cre- 
ated   any    ncticaatle    adverse  effec-s    en   model    fcrecas-s. 

NCGAPS  usc-s  "^he  Araka  wa-S  chu  bert  (1974)  cumulus  paraine- 
-erizazion  ii'-scribed  by  Lord  (1978),  Envircnm  en- al  air 
belofc  the  clcud  -cp  is  entrained  while  air  at  the  clcud  top 
is  detrained.  Ara'^awa  noted  that  large  scale  forcing  gener- 
ates ccnditicral  stability  whereas  cumulus  convection 
destroys  it.  This  concept,  called  •  quasi-equilibriu  ai '  ,  was 
combined  wi":h  the  *  cloud  work'  concept  (Arakawa  and  Schu- 
bert, 1974)  tc  produce  a  cloud  base  mass  flux  unique  to  each 
clcud    type. 

The  GCn  rsdiation  par aireterization  requires  only  fifteen 
percent  of  the  model's  running  time.  The  radiation  transfer 
equation       (Kat  syaiia,  1  97U      and    Sc  hlesinger ,  1976)  uses      bulk 

transmission  functions  fcr  discrete  pressure  layers.  Net 
radiative  flux  2t  ground  level  is  a  function  of  incoming 
solar    and    Icngwave   radiation   and    surface    albedo. 


NCGAPS  is  r'ar.  on  the  CYBER  205  which  has  high  sp^£Q  sca- 
lar ana  vector  pipeline  prccessors.  Presently,  zm  for^oas- 
mo-^.e!  has  a  wall  clock  execution  time  of  thirty  a\inu-e=  p^: 
forecesi  day  en  a  2.4  X  5.0  degree  grid.  Any  significan- 
increass  in  the  model's  performance  would  require  redesign- 
ing the  diata-ic  prccesses  to  include  more  vsctori zation  fo: 
the    en-ire    oicdel. 


19 


III.       SPECTRA  I    DECOMPOSITION    ROUTINE 

Ihe  ccuplex  appaarance  of  a  typical  500  mb  height  wave 
pa-terr.  often  conceals  the  sources  of  forecast  error*  The 
pattern's  C03ipl  =  xity  is  due  to  thr  superposition  of  plane- 
tary, long,  msdiuin  and  synoptic  waves.  Every  pattern, 
regardlsss  of  its  intricacy,  can  be  represented  as  the  sum 
of  a  series  of  pure  harmonic  waves.  A  spectral  decomposi- 
tion routine  will  resclve  a  complex  wave  into  its 
compcnents,  that  is,  the  amplitude  and  phase  angle  of  each 
wave  number  1,  2,  3,  etc.  The  spectral  analysis  used  in 
this  stud  J  was  developed  at  the  Naval  Postgraduate  School  by 
Prof.    R.    I.    wi.lliams. 

The  5CC  mk  wave  pc.ttern  along  U2N  (Fig.  2)  for  18  Febru- 
ary *ill  be  i:pect::ally  decomposed  to  illus-^rate  the 
t  achnique. 

Ic  facilitate  this  illustration  and  the  foregoing 
sequence  of  illustrations.  Fig.  3  has  an  abcissa/or dinat e 
scale    ratic    of    250001crn/60Cm  or   approximately    U2000    to    1. 

Kcte  the  negative  longitudes  are  west  and  the  positive 
lonaitudes    are    east.         A    careful      examination    of    the    heights 
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5G0MB  SFC  HEIGHT 
FEB  18 


90  N 


60  M 


30  H 


SO*  E  90*  E   >20*  E  :50'  E  180*  U  150'  M  120*  M  90*  H   50'  U   30*  u   o* 


T<s»  r   e**  '• 


Figur€    2.       Th€    H^iaht   contcurs   of    the    SOOab   Surface   a-    OOZ 
18    Feb."    Conr.our   Intsivdls    are    120    Meiers. 


HEIGHT  or  500MB  SFC  TOR  42N 


J80.0-150.0  -)20.0  -90.0   -50.3   -30.0    0.0    30.0    50.0    90.0   120.0   ISO. 3   180. ( 

LONGITUDE  I  DEGREES) 


Figure   3.       The    Hiigiit   of    the    SOOmb  Surface    Along    42N    at    OOZ 
18    ret.      The    Nega-^ive    and    Posi-ive    Longitudes    ar 
West    9nd    East    Fespectiv<=ly    of    the    Greenwich 
Meridian.     Boxes   Appear   at    30   Degrae    Intervals. 
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along  42N  suggest  a  relatively  significant  ampli-'id^  of 
wavsnumber  3  as  evidenced  by  the  prominent  ridge  loca-ions 
in  the  vicinity  cf  120W,  0  and  90E.  Also  -he  preponderance 
of  ridges  and  trcughs  in  the  4000km  (50-60  degress)  range 
are  the  3icnatur€S  of  wavecumbars  6  and  7.  Since  a  minimura 
of  six  points  is  required  to  adequately  resolve  a  wave,  the 
maxiiEum  wavenumter  accurately  resclveable  with  a  72-point 
series  {five  degrees  of  Icngitude)  is  12.  k  hcrizontal  r€s- 
cluticn  of  three  degrees  cf  longitude  can  accurately  depict 
wavenumbers  up  tc  20,  Amjlitude  estimates  are  difficult  and 
finer    resoluticn    wcul3   be    speculative    at    best. 

Ihe  nuir€ricai  recults.  of  the  daccmpo  sition  routine  are 
presented  in  lafcie  1  (  although  only  wavenumbers  1-20  can 
be  accurately  resolved,  tire  tables  and  figures  in  this  study 
will  represent  wavenumbers  1  -  28)  .  The  phase  angle  and 
amplitude  (in  rneterb-  about  the  mean  height)  of  the 
wav-E  number    1    through    23    components   are   shown. 

lirst  the  mean  her.ght  was  calculated,  then  th9  amplitude 
and  phase  angle  cf  each  wavenumber  1  -  28.  The  amplitudes 
are  given  as  differ  en  ce  values  abcut  the  m?an  height  and  the 
phase  angles  are  in  r€fer«Ece  to  the  center  of  the  longitu- 
dinal  domain     (0    degrees)  . 
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lAB  LE    I 

Ihe    numsrical  rasuits    cf   the    spactral    decompcsi'^  ion. 

Spectral    decanpositior   of    SOOinb    wave    at  42N.    AraplitudrS   dr=    in 
meters  and    the    phases   are    in    degrees. 


A5IIHMEIIC    MEAN 

WaVE    NUK 

[          AMPLITUDE 

PHASE 

5485.91 

1 

52.04 

-88.94 

2 

66.67 

27.38 

3 

92.17 

-157.89 

4 

24.91 

38.19 

5 

26.40 

126.31 

6 

46.92 

-41.79 

7 

57.17 

129.36 

8 

13.86 

134.82 

9 

13.63 

-0.63 

10 

10.19 

-  115.47 

11 

21.45 

101.96 

12 

10.36 

135.11 

13 

4.26 

-41 .04 

1U 

1.20 

142.10 

15 

6.06 

111.40 

16 

4.96 

125.36 

17 

9.54 

77.27 

18 

4.70 

1.36 

19 

4.54 

64.28 

20 

2.91 

-52.99 

21 

1.13 

-96.49 

22 

1.56 

27.35 

23 

0.46 

-14.74 

24 

1.30 

-146.30 

25 

0.42 

-  114.46 

26 

1.17 

136.61 

27 

1.30 

41.  12 

28 

0.65 

-28.40 

A  graph  cf  wave  number  versus  ampii-ade  (Fig.  4)  fur- her 
clarifies  the  high  relative  amplitudas  of  wavenumbers  3,  6 
and  7  proposed  by  our  iritial  visual  analysis.  Tha  ccmpara- 
-ive    importance    cf    wavenuirber    11     is    also    apparent. 
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SPECTRaL  DECOMF'OSITION  OF  500MB  SfC 


J3.0       23.0       15.0 

WAVE  NUMBER 


Fiaurs   4. 


Ih€    Rclaiive    Anpli'iudss   of    the    Wavenumb^r    1-28 
Compcnsn-s   of    --ne    '500  mb   Surface    Alona   4  2N    at    OOZ 
18    Jeb.    Boxes    Acpeiir    fcr   =  ach    Wavsr.araber    1-28. 


HEIGHT  OF  SOOMB  SFC  flNO  WNl  FOR  42N 


■160.0-150.0 


~i r 

■120.0      -SC.C 


-60.0      -30.0         O.C         30. C        60.0 

LONGITUDE (DEGREES) 


30.0        120.0 


T 1 

J50.0       J80. 


Figure    5.       Th=    Height   of    the    500mb    Surfacs    and    Wavsnumber 
Ccmponents    Along  4  2N    at    OOZ   on    18    F<=t,    The 
Necatlve    a  r.d    Positive    Longitudes   are    West    and 
East    HesoectivEly    of    the    Greenwich    Meridian. 
Bcx€s    are    Positioned    at   30    Degree    Intervals. a) 
Wavenumber    1    Ccjaponent. 
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HEIGHT  or  500MB  SPC  AND  WN2  FOR  42N 


.LONSIIUQCaDJIGBuESI 


5  b)       ttavcnamb^r    2    corapon  en*: 


HEIGHT  OF  SOOnB  SFC  AND  WN3  FOR  42N 


-1 r 

.O-JSC.O    -120-0     -90.C 


-SO.O      -20.0         C.3         30.0 

LONGITUDE (DEGREES  J 


5c)      WGV€n'aaib€r    3    ccEpcnent 


HEIGHT  OF  5C0MB  SFC  AND  WNl-3  FOR  42N 


LONGITUDE  I  DEGREES) 


5d)       Sumffla-ion  cf    wav^namb^r    1-3    compon^nns 
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The  series  ci  graphs  ir  Figs.  5a-e  elucidate  t,he  sigr.if- 
icance  of  the  planetary  ccmpcnent  cf  the  SOOmb  wav?.  All  of 
these  figures  craph  the  height  versus  rhe  longi-ude  cf  -he 
original  wave  (solid  line)  with  various  wavenumber  ccmbina- 
ztonz  (dashed  line)  superposed  for  comparison.  Figs.  5a,  b 
and  c  chart  wavenumbers  1,  2  and  3  respectively,  while  d  is 
the    sumtnaxicn    cf    wavenumbers    1,    2   and    3. 


HEIGHT  or  500M8  SrC  RND  HNl-28  FOR  42N 


-180.0-150. 0    -120.0     -90.0      -60.0      -30.0         0.0         30.0        60.0        SS.O       120.0      150.0       160.C! 

LONGITUDE (DEGREES  J  i 


5e)       Suafiraaticn   of    wavenuicber    1-28    components 


Note  hew  these  three  waves  begin  to  resemble  -^he  origi- 
nal pattern.  It  is  import  ant  to  realize  the  suinraation  of 
the  full  series  cf  harmonic  waves  would  reconstitute  the  500 
mb  wave  excepting  any  high  freguency  variations  (greater 
than    wavenumber    23)     present  in  the  original    wave.         frg.       5e 
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500MB  WN3  PCS!  FOR  42N  NOGRP  0-72 


—   24-0- 


U.1 


t-       48.0- 


72.0 


•ieo.o-:50.3  -120,0  -93.0  -so, a  -3G.o   o.o   3q.o   60.0 

LONGITUDE  I CEGREES) 


90.0       120.0      150.0      la 


U1 


U3 

21 


24.0- 


48.0 


72.0 


500MB  mZ   nNRLYSIS  FDR   42N  NOGflP  0-72 


-1BO.0-I5O.3  -120.0  -SO.O 


-60.0   -30.0    0.0    30. C    60.0 

LONGITUDE [DEGREES) 


30. 


150.0     la 


Figure   6. 


Hovircil*^-    Diaqrams    of    tha    Wavenumber    3    Compon<?r- 
Af.aiyses    ara    fcr  OOZ    en    16-21    Fbd    from    NOGAPS. 


0 

liirervais    are    3C   Ma-srs 
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rspres^nes  the  suid  of  waverumb^rs  1  through  28  and,  qt.  -his 
scale,    appears    tc   coincide    wir.h    th-3    zss"    wave. 

Ihis  analysis  approach  has  prov-id  zo  fce  a  useful 
research  tccl  when  applied  t.o  fcracasr  verifica-ion  (Baum- 
hefner  and  Dcwney  ,1  978;  Scm^r  vii  le,  1  980;i  .  Specxral  decompc- 
si-icn  of  a  forecast  and  iz s  verif icition  analysis  allow 
direc-  cotppariscr.  of  rhe  aipiitudrr  and  phase  speed  of  each 
wavenumber.  The  dimension  of  time  fur"her  sireng-hens  "his 
analysis    "technique      via    a    Hovaioller   diagram.  These    trough 

and  ridge  Icngi-uca-i;  ime  plots  along  a  fixed  latitude  depict 
the  time  evcluticn  of  the  forecast  variables.  ?ig.  6  ccni- 
pares  the  anclitude  and  prase  speed  of  successive  24,  U8  and 
72  hour  forecasts  for  wavenumber  3  with  th=  observed  atmos- 
pheric  behavior. 

In  the  next  chapter  this  and  other  cases  will  be  ana- 
lyzed in  this  Banner.  The  hamorics  of  each  forecast  are 
grouped    into      planetary  (1 -3)  ,       Icng    (5-7)  and    medium  (8- 12) 

waves  befcre  comparison  to  the  verification  analysis  via 
Hovmcller    diagrams.  Headily    discernable    errors      in    ampli- 

tude, phase  speed,  develcfiient  of  baroclinic  systems,  damp- 
ing cf  smaller  scale  features  and  planetary  wave  struc"^ure 
are  just  a  few  of  the  advantages  cf  verification  using  spec- 
tral   comocnents. 
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IV.       SPECTRAL    VEJIFICATION    CkSE    STUDIES 

In  rhis  chapter  the  spectral  verification  presented  in 
chapter  3  \iill  te  applied  to  thrs?-  NCGAP5  forecas-s,  one  for 
three  days  and  two  extending  for  five  dc.ys.  All  forecasts 
are  lade  from  OOZ  data  anc  are  verified  c.t  24,  U8,  72,  96 
and    120    houi    lerks  of   the    forecast. 


4.1      Case    1    -    Fetruary    18    tc    21,     1982: 

This    period    was    selected   on      the    bas:.s    of   its    crcncuncsd 


planetary    and    long  wave    activi-y.        ''"his    prediction   was    m; 


^.- 


with  an  earlier  »coarse  lesoluticn'  (4  degrees  lat.  by  5 
degrees  long.)  version  of  the  NOGAPS  f:;recast  model.  Tha 
spectral  decomposition  was  applied  to  features  along  42N 
averaged    k€tw€sn    38N    and    46N. 

The  initial  SOOmb  pattern  for  case  1  was  discussed  in 
Chapter  3  and  is  shown  in  Fig.  2.  Prominent  features  alone 
42N  include  planetary  wave  troughs  at  160W  and  45w  and 
strong  ridges  at  170E,  1  lOW  and  15E,  suggesting  a  strong 
wavenumber      3    ccirponent      aiplitude.  The    preponderance      of 

ridges  and    troughs   with    wavelengths    on  the    order   of   50    to    60 
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deqr€6s  are  -he  signatures  of  wavrnambars  6  and  7.  B^call 
zhat  spactial  decomposition  of  the  heights  along  42N  con- 
firmed the  relative  dominanc'r  of  wavenumber  3,  7  ?.nd  11 
components    (Fig.    4)  . 

The  72-hcur  forecast  and  analysis  for  the  SOOmb  surface 
on  21  Feb  (Figs.  7a  and  b)  appear  quits  sinilar  under  a  cur- 
sory   inspection. 

The  trough  at  60-70  W  has  been  under  fcracc.s-  by  the  model, 
but  overall  the  forecast  seems  to  be  reasonably  successful. 
Note  the  excellent  forecast  of  the  two  troughs  in  the 
Pacific  Ocean  (180W  and  1U0W)  and  the  5-trong  ridge  over 
Europe. 

iihen  the  spectral  deccraposi ticn  is  applied  to  the  72 
hour  forecast  and  its  verification  (Fig.  7c)  significant 
amplitude  differences  are  present.  opec:.f  ically ,  in  the 
planetary  wave  group  (wavenuabers  1-3)  tha  observed  atnios- 
phere  has  maiitained ,  although  to  a  lesser  degree,  the 
dominance  of  wavenumber  3  whereas  the  model  has  shif-^ed  the 
peak  amplitude  away  from  vfavenamber  3  and  into  wavenumber  1. 
In  the  long  wave  group  (wavenumbers  5,6  7)  the  dcminanc=  of 
wavenumbers  6  and  7  has  shifted  to  the  sole  peak  in 
wavenumber    7    in    both    the    iicdel  and   the   analysis. 
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a) 

90*  M 

SO*  H 
30*  N 


0* 


500MB  SFC  HEIGHT 
NOGRPS  72  HR  FOR  FEB  21 


60*  E  90*  E   120*  Z     ISO*  E  190*  «  ISO*  U    J20*  H  90*  W   60*  M   30*  W   0*     30*  E   60*  E 


b) 


90  N 


SO*  H 


30*  H 


500MB  SFC  HEIGHT 
FEB  21 


60*  E  90*  E   120*  E  ISO*  E  180*  H  ISO*  W  120*  M  30*  W   60*  W   30*  M   0*     30*  E   60*  C 


FiQur?    7.       SOGmt    Heights    and    Wav<^niimber  Ccoipor-entf      for    OOZ 
'  21    Pet.      a)    Forecas-     ana    o)     Anaxysis    w_-h   Cor.  .cu. 

Intervals    cf    120   ^.eteis.    c)  Wavenumber    Ccmpon3n- 
Ampiitudes    in    iiie-csrs. 
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c) 
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SPECTRAL  DECOMPOSITION  OF  500M8  SfC 

RNflLYSISISXID)  fiM3  72HR  PORECflSTlORSHED; 

21  FEB  AT  42N 
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WRVE  NUMBER 

7c)    Wsver.umber   coaiponsr.t    aaplitudas. 


NOGAES  accura-€ly  for-cast  nh^  wavanamDer  a.iipli'-.ud  =r-  in  -his 
group.  Th€  Eediuai  wave  group  (wavenumber  8-12  ccn?por-€  r.  ts) 
initially  indicates  a  prafsr^ncs  fcr  wavenuraber  11,  tut 
after  7  2  hcurs  -.he  significance  has  shift  ad  zo  "he  lower 
wavenunibers  8  and  9.  This  shift  is  slightly  mor'?  proncunced 
in  the  nodel  resulting  in  erroneously  high  forecas-:  a.iipli- 
tudes    for    waverumfcers   3    aid   9. 


Iha    Hovmoller   diagrams    for   the   planetary    wave    ccmpcnents 
of    the    forecast    and   the    analysis       (Figs.       8a    and    b) 
illustrate    these    differences. 


,- 1  c  = - 


J- 7 
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a) 


PW  FCST  FOR  38-4oN  NOGRP  0-72 
FEB  18  -  FEB  21 


300.0  270.0      240.0      210.0 


180.0       ISO.O      120.0       90.0        6C  .0 

LONSlTUDEt DEGREES) 


30.0    0.0   -30.0   -50.0 


b) 


PW  RNRLYSIS  FOR  38-45N 
FEB  18  -  FEB  21 


300.0  270.0      210.0      210.0      180.0       150.0      120.0       33. G        60.0        30.0         0.0 
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-30.0      -50.0 


^iqur€   8.       Planetary    Wava    Hcvmoller    Diagrams    for    18-21    ?<=-h 
a)    rcr€cast    and   b)     .\nalv3i3    with    Contcur 
Ir.rsrvais    of    5C    Me-.-rs.' 


The    iDcdel      has    accurately    positioned      ridges    at    120W      and   0, 
bat    it    hao    undsres-imated    the    amplitude      of    the    former    by    50 
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meters  and  cversstiraa ted  the  latter  by  50  meters.  The  -wc 
planetary  wave  trcughs  criginaring  at  180  W  and  ilOW  ar  ^  weak- 
ened too  rapidly  in  the  forecas"^  and  becDme  distorted  by  the 
spurious    fcrms-icn      of   a      deep  trough      at    9CE.  The    trough 

originating  at  30E  also  weakened  too  rapidly  -o  -h-?  point 
where  iz  is  higher  rhan  the  mean  heigh-  and  is  only  a  trough 
in  ths  sense  -ha-  i-  is  bstween  two  ridias.  The  ct.her  prom- 
inent feature  is  the  forecas-  trough  iz  ?0E  which  has  no 
counterpart      in    the      observations.  The      iieemingly      randcni 

changes  in  trough  anc  ridge  intensity  for:i  a  pattern  wh^n 
the    differences    m   component    ampiitud'is   are    considered. 

The  NCGAFS  72-hour  forecast,  shows  the  -^avenumber  1  sig- 
nature emerging  with  general  ridging  from  150W  to  UOE  and 
ganeiai  trcughing  from  UOE  to  ^50^  while  the  analysis 
depicts  mors  cf  a  wavsnumter  2  and  3  signature.  Figs.  9a 
and  t  compari  the  forecast  and  cbservel  wavenumbers  1  and  3 
inai^idua  lly. 

NOGAcS  forecasts  wav^nuinter  1  to  be  70S  more  intense  than 
the  cbservation  while  underforecas- ing  wavenamber  3  by  55-^. 
Despite  the  presence  of  wavenuaibers  2  and  3,  the  planetary 
wave  char2C-:€ristics  in  the  72  hour  NOGAPS  forecas-  are 
largely      determined         by      the         overwhelming      dominancr        of 
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a) 


flNRLYSIS  WNl  COMPONENT  ("SQLTDT 
FORlCAST  WNl  COMPONENT  (OaSHEO) 
72  HOURS   (21   TEBJ 


30C.0  270.0      240.0      210.0 


180.0       350.0      123.0       30.0        5J.0 
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■30.0   -50.0! 
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b) 


BNRLYSIS  WN3  COMPONENT  150L1DI 

FGRECflST  WN3  COMPONENT  lOBSHED) 

72   HOURS  121  FEB) 


300.0  270.0      240.0      220.0 


180.3      150.0       .20.0       90.0 

LONG I TUOE I  DEGREES} 


Figure   9.       iia  v  snuaber    1    ard   3    Co  mpcn^rxts    for    Case    1.       a) 

Waver.umber  1  Ccirponent  for  Analysis  (Solid)  and 
Foiecas-  (Cashed)  for  OOZ  21  rso.  b)  Same  as  a 
excep-    for   Waver.umbsr    3  Component. 


waver.umbsr  1.  Zrronacusly  lew  heights,  troughs  too  strong 
and  ridges  toe  weak,  exist  in  the  region  of  60Z  to  120i^  and 
erroneously    high    heights,        troughs    too   weak    and      ridges    tec 
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Swrcng,  exis-  in  the  regicn  of  1 2QW  to  60E.  Synoptic  sxam- 
pl€S  cf  the  effect  of  the  raduction  of  the  (*av  er.uaber  3 
compcnent  ere  evident  in  the  ridge  a-  120W  and  -.h-  trough  at 
60W.  From  150W  to  1 10W  the  model  under  forecasts  the  ridge 
in    this    area      and  smooths    cut    its    ainpl:.tude.  From    110W    to 

60W  the  model  over  forecasts  the  height  of  the  aastern  side 
cf    the      ridga    and      smooths    cut      the    trough      at    60W.  These 

amplitude  errors  are  also  evident  in  other  regions  of  the  72 
hour    forecast    (Figs.    7a    aid  b)  . 

Ihe  Icng  «ave  group  verifies  the  wav?number  5,6  and  7 
compcnents,  Th=  general  agreement  betw^^en  tiie  individual 
wavenumber  aaplitudes  in  this  grcup  (Fig.  7c)  through  th^ 
72-hcur  period  ccnseguently  produces  a  ge!:eral  agreement  in 
the  airplitudes  of  their  respective  Hovmoiler  diagrams  (Figs. 
10a    and   b)  . 

This  diagram  shows  seven  troughs  an i  ridges  at  72  hears  sig- 
nifying the  dcninance  of  viavenumber  7,  in  this  grouping  a 
regicn  of  maximum  amplitude  (contour  interi'als  are  15m)  is 
eviden-  between  180W  and  SOW  in  botn  tha  model  and  the 
observation.  There  is  strcng  agreement  between  the  forecast 
and  observed  long  waves  in  the  Hovmoiler  diagrams  shewing 
NOGAFS    Derfcrmina  well   at    this   scale. 
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LW  rCST  FOR  38-46N  NQGRP  0-72 
'   FEB  18  -  FEB  21 
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Figure  10.  Lor.q  Wave  Hovmcliar  Ciagraras  for  13-21  Feb.  s) 
Fcrecasr  and  t)  Analysis  with  Coniour  Intervals 
of    30   Me-ers. 


Ihe      in-eiac-ion    of       v»avenum  ters      8      -hrough    12      in      the 
medium   wave    grcuf  creates    a      complex    pattern    in    the   analysis 
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and  forecast  trough-ridge  diagrams.  Tha  Hovmoller  diagrams 
for  the  fcreccst  and  observation  (Figs.  11a  and  b)  both  show 
eleven  waves  initially  and  nine  waves  af-er  72  hours. 
The  airplitudes  (contour  intervals  are  10  inerers)  ar.  d  phase 
speeds  can  s-ill  he  compared,  but  zh3  errors  ar-s  mor-r  diffi- 
cult to  trace.  The  model  shows  an  erroneous  less  of  ampli- 
tude, particularly  evident  ai  US  hours.  Hcw^vrr,  a-  72 
hours  the  increased  activi-y  in  the  region  frofi?  140E  to  120W 
is  properly  forecast.  The  phase  speed  is  aporoxima- ely  10 
degrees   -co      fast   in      -he    region      fr^m   00 E      tc    120Z.  Con- 

versely, the  model  has  done  extremely  well  in  f cr'^cast-ing 
the  speed  of  the  wave  features  between  li+OE  and  I50W.  Over- 
all, The  mccel  exhibits  reasonable  s\icc'=:ss  for  both  the 
amplitudes  and  phase  speeds  cf  the  medium  waves  a-  72  hours. 
I-  is  import anr  tc  no-e  the  amplitude  of  -he  aedium  rfavenum- 
bers  is  often  smaller  than  the  planetary  wavenumbers  by  as 
much  as  an  order  of  magnitude.  Also  as  wavenunbers  increase 
to  the  point  where  they  cannot  be  accurately  represented  by 
the  model's  horizontal  resolution,  the  spectral  decomposi- 
tion   is    less    accurate  as    fcell. 

The      results      of      Case      1      indicate      NOGAFS      erroneously 
shifted   the    planetary      wavenumber   dominance    from    wav=      3    and 
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Figure   1 1  . 


Medium  'rfavs  Hcvmciiar  Diagrams  for  18-21 
a)  F  CIS  cast  ard  b)  Analysis  wi-h  Conto-jr 
I.ixervals   of    15   M'srsrs. 


-30.0      -60.0 

Feb. 


into    have    1.         Du?  -o   a    difference   in    wavenumber   signa-.ures, 
accurate    planetary   wave    phase      speed   comparisons    beiween    the 
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model  and  tha  observation  were  difficult.  Both  the  ampli- 
tudes and  phase  speeds  cf  zb.<=  long  waves  were  iccaritsly 
predic-ed  ty  the  model.  NOGAPS  shewed  some  less  cf  ampli- 
tude in  the  ledium  waves  and  shifted  peak  amplitude  from 
wavenumber    11      to  8    and    9.  The   phase  speed   errois      of   the 

mediuir   waves    as   a   group    were   small. 


4.2      Case    2    -    October   06    tc    11,     1982: 

This  pericc  provides  a  good  example  of  the  fcrTiaticr.  cf 
an  intense  ridge  (blocking  high)  over  the  western  UcS.  and 
was  discussed  in  the  NOAfi,  Western  Regional  Attachmen*  of 
Octot^r ,  1 S82.  For  this  case  NOGAPS  horizontal  r;-3clii*ior  is 
2.4  degrees  of  latitude  ty  3.  degrees  of  longitud?  znterpc- 
lated  onto  a  2.5  cy  2.5  degree  spherical  grid.  The  ^pec-ral 
decompcsit  ic  r.  is  applied  tc  the  height  of  the  SOOnih  sur!:ace 
along  45N  (avaraged  tetween  the  heights  given  at  42. 5N  and 
47  .5N) . 

The  50Ciiit  heights,  depicted  from  30N  to  60Ny  for  06  Cct 
in  Fig.  12a  shows  a  fairly  zonal  pattern  fr^m  the  east  coast 
of    Asia   to    th€    Easter r.   Atlantic. 

Ncte  the  ordinate/abscissa  ratio  is  mcr^  than  one.  This 
distcrtion    produces    a    visual    magnification    cf    the    amplitudes 
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a) 


500MB  SFC  HEIGHT 

OCT  06 


60*  E      90*  E        120*  E       ISC*  E       :B0*  W       550*  H       120*  W       90*  W        60*  W        3C*  H         0*  30*  E        £ 


b) 


SPECTRRL  OECOMPOSiTION  OF  500M8  SFC 
flNfiLYSlS  08  OCT  AT  45N 


Fiaure    12. 


WfiVC  NUMBER 


5CCint  Heiahts  and  Wavsr.umber  Cornpori-3nr: 
Cc--  a)  SOOmt  Heights  wi-'n  Con*  our  In- 
120  Meters,  b )  Wavenu  aifa?r  Ccmponen~  Amp; 
Me-ers. 


*cr    06 
rrvils    of 
.tudes    :.n 


of  each  wave.  This  scale  was  chosen  for  clari-y  and  is 
maintained  -hrcughout  this  and  the  remaining  cases.  Promi- 
nent   ridaes    are      loca-ed    at   20W    and   30E  and      -he    flow    scli-s 


ai 


at  10W.  Spectral  deccmposition  (Fig.  12b)  shews  wav^uanifc^r 
2,  5  and  7  compcnents  ar£  dcminar.r  ini-ially.  Th?  zonal 
character  cf  this  case  is  shown  by  the  relatively  even  dis- 
tribution   cf    anplitudeo    ancng    wavenumbers    1    through   8. 

S€ventj-t\iO  hours  later  (Fig.  13a)  a  deep  -.rough  has 
developed  ever  the  central  U.S.  and  prominen-:  ridges  have 
rapidly  develcped  over  the  western  and  eas-::arn  U.S.  ccusts. 
The  model  (Fig.  13b)  does  not  adequataly  forecast  the 
strength  cf  these  features.  A  ccmparison  of  the  sp-.c-'ral 
decompositio r  (Fig.  13c)  cf  the  analysis  and  the  forecast  at 
72  hcurs  shows  a  marked  increase  in  the  dominance  of  the 
wavenuaber  2  and  7  ccmponents  of  the  atmosphere  whil?  the 
mcdel  incorrectly  smoothed  the  amplitude  of  wavenumber  2  but 
correctly    mairtained    -^he    relative    strength    of    wavenuraber    7. 

At  the  120  hour  mark  (11  Oct)  of  the  forecast  period 
both  the  eastern  and  western  U.S.  ridges  have  ccr.tir.uea  tc 
strengthen  (Fig.  14a) .  Ihe  model  shews  an  ability  to  prop- 
erly forecast  the  phase  speed  of  the  wave  features  such  as 
the  troughs  at  170W  and  40W  and  the  ridge  at  160E  but  con- 
tinues to  underestimate  the  strength  of  the  ridges  (Fig. 
I4b)  over  the  U.S.  Spectral  decomposition  of  the  observa- 
tion   shows    the   largest    amplitudes    are    in    wavenumbers    2,    3,    U 
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a) 


500MB  SFC  HEIGHT 
NOGRPS  72  HR  TOR  OCT  09 
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C) 


SPECTRAL  OECOrtPOSITION  Of  500MB  SfC 
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09  OCT  AT  <15N 
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WAVE  NUMBER 
Fiqure   13.      Same    as    Fig. 7    -except    for    09  Oc*. 
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and  5  all  having  similar  valaas  (Fig.  14c).  M?5ar.wh.il=  -h« 
model  has  shifted  the  dominance  to  wavenamber  'i  whii^  kz  3p- 
ing  the  amplitude  of  wavenumbers  5  and  greater  reiativ^-ly 
small.  The  largest  amplitude  changes  in  the  mo^^el  during 
the  forecast  period  occurred  in  wavenumbers  2,  3  and  U  whiis 
the  atmosphere  additionally  showed  large  amplitude  changes 
in    wavenumter    5. 

Ihe   sicnificancs    cf    the  amplitude    errcrs    is    readily    di=- 
carnafcle    in      the    Hovmcller     plots.         Initially      bo-.h   of      -hi 
planetary      wave    troagh-ridge      diagrams      (Figs.         15a    and      b) 
exhibit      the        characteristics      of      a      large         wavenumber      2 
ccmpcnent . 

After  "72  hours  the  observation  has  greatly  increased  tb~ 
strength  cf  the  wavenumber  2  signature  while  the  forecast 
decreased  it  <Fig.  16a).  The  model  has  successfully  fore- 
cast the  strength  of  the  trough  at  180W  (-135in),  but  the 
strength  of  the  other  trough  and  both  ridges  are  severely 
underestimated.  The  combined  dominance  of  wavenumbers  2  and 
3  in  the  NCGAES  forecast  has  introduced  a  spurious  ridge  at 
OW    and   trough    at    50  2. 

Ey  120  ncurs  the  atmosphere  exhibits  maximum  aapli-ude 
in    the    planetary    wavenumbers    2   and  3     (Fig.    14c) .       This    shift 


44 


a) 


500MB  SFC  HEIGHT 
NOGRPS  120  HR  FOR  OCT  11 
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SPECTRAL  DECOnPOSITION  OF  SOOMB  SPC 

flNflLlfSIS(  SOLID}   AND   120HR  FORECAST  (DASHED) 

11  OCT  AT  45N 
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Figure    14.      Sam*;    as    Fig. 7    except    fcr    11   Oct. 
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Figure    15.      Same    as    Fig.     6   except   for   06-11    Oct. 

has  net  affected  the  amplitudes  cf  the  planetary  wave  fea- 
tures as  much  as  it  has  sharply  tightened  the  spacing 
between      the      trcugnes      ard      ridges         producing      more      cf      a 
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wave  numb -5  1  3  signature.  The  model  has  forecast  -he  aacii- 
tude  of  waver-umter  1  accurately  while  under-?stiina- ing 
waverumbcr    3    by    30%    (Fig.     15b). 

Ihus  ths  model  has  cnce  again  failed  to  predict  th« 
shift  and  errcnecusly  smocthed  the  highest  amplitudes.  The 
model's  planetary  wave  Hcvmoller  diagram  shows  the  weak'=^r 
amplitudes  at  120  hours  and,  significantly,  thr  rssuliant 
smoother  pattern.  Specifically  the  ridges  at  105E  and  110W 
are  tco  low  and  the  trough  at  '45W  is  tec  high.  Synoptically 
th=se  errors  are  evident  ir.  Fig.  lUb  as  the  failure  of  th  = 
model  to  adaquately  forecast  the  strength  of  -he  ridge 
building  on  *he  west  coast  of  the  U.S.  and  the  trough  at 
40H.  The  general  placement  of  comparable  features  is  gccd 
especially    in    the  Pacific    region    near   the    dateline. 

The  long  wave  component  Hovmollar  diagrams  (Figs.  17a 
and  t)  shew  NCGAPS  successfully  forecasting  the  long  wave 
features.  Overall  the  forecast  has  accurately  predicted  -he 
general  regions  of  high  long  wave  amplitudes,  120W  to  90E, 
and  lew  long  wave  amplitudes,  at  902  to  120'.*.  Observed  long 
wave  amplitudes  are  more  intense  than  the  forecast  after  the 
72  hour  mark.  Comparison  of  model  and  observed  phase  speeds 
shows   NOGAES    to    be   fast,       particularly   m    the    high   amplitude 
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Figure    16. 


rtav  =  r.uinber   2,    3   ani    5    ComDcnsn- s    for   Case    2. 
Sai€    as    Fig.     9   excspt   a)     Wave.iumbsr    2    for    09 
Oct,    b)     Wavenumber    3    for    11    Oct    ani    c) 
Wavenumber   5    for    11    Oct. 
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wavss  between  120W  and  30W.  The  shift  of  -he  atraosph&re  -c 
wavsEuraber  6  signature  xitnin  the  group  by  72  hc'irs  was 
properly  forecast  by  the  modal  and  wavenumber  5,  6,  ard  7 
components   hav5   less    than    10%  ampli-cude  error    (Fig.    13c). 

Ey  12C  hours  the  atncsphere  shows  maximum  long  wave 
amplitude  in  the  wavenumber  5  component  (Fig.  lUc)  and  the 
Hovmcller  clot  (Fig.  17a).  In  contrast  the  model  (Fig.  17b) 
exhitits  snail  long  wave  amplitudes  and  underestimates 
wavenumber  5  by  more  than  50%  (Fig.  16c)  and  wavenumber  6  by 
30%  while  overestimating  wavenumber  7  by  60%.  The  resul*  of 
thes«  errors  is  a  general  smoothing  of  all  ridges  snd 
troughs.  Also  the  region  of  high  amplitude,  120W  to  90E, 
has  exhibitid  slew  phase  speed  error  while  the  region  cf  low 
amplitude,    90r   to    120W,    shewed   rapid    phase    speed. 

I he  medium  wave  components  (Figs.  18a  and  b)  begin  the 
forecast  period  with  the  signature  of  wavenumber  8.  By  72 
hours  no  clear  viavenumber  preference  is  present  in  either 
the  forecast  or  the  atmosphere.  The  same  is  -^rue  after  120 
hours.  The  initial  significant  amplitudes  of  these  waves 
ar€  lost  ir.  the  atmosphere  after  24  hours.  A  few,  such  as 
the  trough  originating  at  lUOW  and  the  ridge  at  120W,  main- 
tain   their    amplitude    throughout    most    cf   the    five    day    period. 
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Figure    17.       Sania    as    Fig.     10   except    for    06-11    Oct. 


Az  120  hcurs  new  activity  appears  in  the  region  from 
130W  to  SOW.  The  model  dampens  all  the  wave  fea-ures  afrer 
2U    hcurs    «hich   is  errcneous  in      the    case   of   the   trcugh-ridge 
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Figur?    18.       ScIds    as    Fig.     11    except   for   36-11    Oct. 

sys-€!D  over  the  western  U.S.  The  modal  doss  capture  the  new 
activity  at  the  five  day  mark.  An  Dverali  ccmpariscn  sug- 
gests  a      phase    speed      errcr  of   5-10      degrees    too      fast    where 
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direcx  comparisons  can  be  made.  Synopt.  icaily  this  differ- 
ence is  nc^  visible  en  the  scale  of  Frgs.  lUa  and  b  tut 
would  shcvi  up  as  a  slightly  rapid  aiovemen-  of  fea-ures  with 
a  wavelength  cf  30  to  U5  degrees  over  zhe  120  hour  forecas- 
ter icd. 

The  results  of  Case  2  indicate  NOGAPS  erroneously 
smoothed  a  shift  in  the  planetary  wave  act-vi-y  at  72  hours 
and    again    ct    120    hours.  Planetary    wave    phase    speed    errors 

were  small  fcr  comparable  features  and  negligible  in  the 
Pacific  region.  Although  long  wave  forecast  amplitudes  were 
smoothed,  MCGAPS  correctly  produced  two  distinct  regions  cf 
contrasting  anflitudes.  Phase  speeds  were  slightly  too  fast 
in  the  high  amplitude  region  but  satisfactory  in  the  low 
amplitude  region.  The  airplitudes  cf  the  medium  wave  fea- 
tures were  well  forecast,  but  shewed  a  fast  bias  in  wave 
speed    and    several   Key    features    were    lost   m    the    forecast. 


4.3      Case    3    -    Cctcber   22    tc   27,     1982: 

Curing   this      forecast    period      the    500mb      surface    changed 
from   a   predominantly    long    wave   flow      pattern   to    a    more    zonal 


situation.  In    particular      the    region      from       150 


lUOW 


underwent    the      brea)cdcwn    cf      a   bloclcing      ridge.         Again      the 
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NOGAES  hcrizcn-^.al  resolution  is  2.4  degrees  la-.  by  3 
degrees  long.  interpolated  onto  a  2.5  by  2,5  degrse  spheri- 
cal grid.  The  spectral  decomposition  is  applied  to  the 
heights    alorc    U5N   averaged   between  42. 5N    and    U7.5N. 

Ihs  500inb  heights  from  30N  to  60N  for  22  Oct  (Fig.  19a) 
includes  prcmirer.t  troughs  at  60E,  125E  and  140W  with  lesssr 
troughs    at    IS'a   and    15W. 

Noteatle  ridgss  exist  at  90E,  170E,  110W,  45W  and  202.  The 
wavenumber  5  dcminance  is  clearly  present  in  Fig.  19b  aicnc 
with  strong  aiplitude  compcnents  of  wavenumbers  3  and  4  due 
to    their    particularly   3t.rcng   trough-ridge    systems. 

After  72  hours  amplitude  differences  are  already  appar- 
ent tetween  •'•he  500mb  patterns  of  the  model  and  the  analysis 
(Figs.  20a  and  b) .  The  model  has  underfo recast  the  strength 
of  the  ridges  at  90E  and  25E  and  the  troughs  at  70E,  14GE, 
140W    and    ICE.  The    only    obvious    NOGAES      ovarforecast    is    in 

the    lidge    at    100^. 

The  individual  wavenunber  components  (Fig.  20c)  a-  72 
hours  show  the  dominance  remains  with  4  and  5.  Larga  ampli- 
tude increases  have  occurred  in  wavenumbers  1,  7,  9  and  10 
and  large  cecreases  at  2  and  3,  The  model  has  accurately 
predicted         most  of        the        changes  but         ccnsiste  r.tly 
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WftVE  NUMBER 
Figure    19.       Same    as    Fig.     12  excep-    for   22    Oc- 

ur.  der€stiira-  =  s   tht    strength   of   wavenumbsr  s  6    through    10.       In 

tarms   of    this      study's   wave   groupings    the      amplitudes    of    th.= 

planetary    wavas    are    forecast   very      well,      the   long   waves    are 

too    weak    ty    10      -c   207c   and   the   medium    waves      are   too    weak    ty 

30     tC     UO'':. 
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Figure    20.      Sams    as    Fig.     7   except   for    25    Oc- 
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at  th€  120  hour  mark  both  pha  s-r  spead  and  amplitude 
errors  are  svidsnt  in  the  SOOmb  patterns  (Figs.  21a  a^.d  b)  . 
Phase  speeds  are  consista  rtly  tec  rapid.  Examplis  include 
the  trough  at  75E  (10  degrees  of  longitude  too  fast),  the 
ridge  at  180W  (15  degrees  too  fast)  and  the  trough  at  uOW 
(25  degress  tec  fast).  The  amplitudes  of  th^*  trough  at  160S 
and  the  ridge  at  180W  are  too  weak  and  the  ridge  at  95W  is 
too  strong.  A  ridge  has  formed  in  the  ragion  from  160W  to 
120W  which  has  aided  in  the  weakening  of  the  blocking  ridge 
at  1S0W.  The  model  has  missed  this  important  feature 
entirely. 

The  individual  wavenuaber  components  (Fig.  21c)  at  +he 
120  hour  mark  shew  a  general  decrease  in  all  amplitudes  wi*h 
especially  large  percentage  drcps  m  wavenumbers  3,  U,  5,  6 
and  7.  The  relative  dcainance  of  planetary  waves  is  much 
smaller  new,  tut  there  exists  a  distinct  peak  at  wavenumher 
U.  The  mcdel's  wavenumfcer  components  also  show  a  general 
decrease  in  amplitudes  with  significant  errors  occurring  in 
3,  H,  8,  S  and  12  while  the  dominant  wavBnumber  is 
wavenumber  3.  In  terms  cf  this  study's  wave  groupings  the 
amplitudes  cf  the  planetary  waves  are  too  strong  by  20  to 
30%,  the  long  waves  are  wsll  for<=cast  and  the  medium  waves 
are    tco    weak    ty    3  0  to    H0%. 
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Figure   21.      Sam?    as    Fig.     1   except    for    27    Oc- 
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1h<5  Hcvnicll  =  r  diagrams  for  th=  plar.etary  wave  corr.poner.t  s 
of.  -the  forecast  and  the  aralysis  (Figs.  2  2a  and  b)  both  show 
a  shift  from  ths  sole  dominance  of  wavanumbar  3  zc  a  shared 
dorairance    with    wavenumber    1    by  7  2    hoars. 

Ths  forecast  pha3€  speed  error  of  -he  model  is  r.egligibl=  a- 
this    point.  The    wavenuirber      3    f=^atur=s       are    forecast      too 

s'lror.g  as  a  result  cf  the  20%  ovarf orecasi:  by  the  model. 
Examples  include  the  troughs  a-  1102  and  150H  and  the  ridges 
at  8CS  and  30E.  A-  120  hours  the  model  shifts  i-s  maximuni 
activity  tc  wavenumber  3  while  the  atmosphere  shews  a  dis- 
tinct wavenumber  1  signature.  Now  the  NOGAPS  wavenumber  3 
features,  the  troughs  at  120E,  SOW  and  35E  and  the  ridges  at 
85W  and  35E,  are  much  too  strong.  Fig.  23  graphically  pres- 
ents the  lOC^  cverforecast  of  the  wavenumber  3  componen-  by 
th=    ircdel. 

The  pcsiticns  of  the  ridges  at  85W  and  35H  are  excellent, 
but  the  trough  at  120E  has  been  erroneously  retrograded  15 
degrees.  Another  interesting  aspect  of  this  case  is  the 
retrogression  cf  the  planetary  wave  which  originally  was  at 
150W.      NOGAES    handles   this   phase    shift   very    well. 

Within  the  Icng  wave  group  the  signature  of  wavenumber  5 
persists    through      72    hours      in    both      the    model      and   analysis 
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Figure    22.      Same    as    Fig,     8   excspt   for    22-27    Oct 


(Figs.  242  anc  t)  .  A  general  ccmparison  cf  obserV'=:d  and 
forecast  Icng  wavss  shows  good  agreeroint  on  this  seals, 
particularly    in    phase  speeds.         The    model    anderf  or  Tcas-s    -h^r 
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'igu]:€    23 


Wave  number    3    comoonan-    for   Cas9    3 
9    €xcsp-    Wav<=  number    3    for    27   Oc-. 


Same    as    Fig. 


~':r  =  !:g-h      cf    wavsr.uinbars      6  and      7      producing    sligh-.ly      w^ak 
s..Tipl:.tud5s    in   th€   region    frcm    30S    -o    120S. 

However,  by  120  hears  both  observed  and  forscasz  TDa-.- 
"ern.'S  show  a  duel  dominance  in  wave  numbers  5  and  6  with  the 
jiodal  amplitude  a  li-tle  too  weak  in  wavenumber  5  and  a  li-- 
-le  too  strong  in  wavenumter  6.  In  -he  iasz  2U  hours  of  the 
forecast  c€ricc  the  atmosphere  introduced  a  six-h  trough  and 
ridge  by  a  general  retrograding  of  certain  existing  features 
(180W  to  UCE)  by  5  to  1C  degrees  of  longitude.  The  model 
introduced  a  sixth  trough  and  ridge  by  a  rapid  eastward 
movement  of  certain  features  (25W  to  1U0E)  .  This  incorrect 
me-^hcd  of  hcncling  the  dcminance  shift  resulted  in  phase 
speed  errors    in   ■^he    range    of   5  to    15    degrees    too   fast. 
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Figure   24.       Same    as    Fig.     10   excsp-    for    22-27   Oc- 


In  -he  madiurn  wave  group  (Figs.  25a  and  b)  -he  model 
correc-ly  began  with  -he  signature  of  wa  ver.uniber  9,  .uain- 
-ained      Lz    through      72   hours      then   shifted      to    wavenumber      3 
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Figure    25.      Sam^    as    Pig.     11    except  for   22-27    Oct 


doaiir.ar.es  ty  120  hours.  Both  diagrams  show  th^ir  largest 
amplitudes  in  ths  region  from  30 E  to  120W.  Overall  the  mod- 
el's   amplitudes      are    too      weak  as     a    result      of    consistently 
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und«=iicr^cas-ir.g      the   strength      of   the      medium    waves.  The 

phase  speeds  are  too  fast  by  apprcximately  10  degrees  in  the 
regicn   frcni    2CE   tc    16 OE. 

The  results  of  case  3  indicate  NOGAPS  erroneously 
shifted  the  planetary  wavenumber  dominance  into  wav<=numter 
3,         The    error    «as   large      enough    to    significantly    effect    the 

NCGAPS  plane- 
tary ^ave  phase  speed  was  excellert  at  72  hours  and  still 
good    by       120    hours.  The   forecast      acplitudes    of      -he    long 

waves  were  only  slightly  sttcothed  and  long  wave  ph'.se  speeds 
were  exc  =  ll'T:nt  at  72  hours  but  too  fast  by  120  hours.  The 
mediuic  wave  airtlitudes  were  smoothed  and  phase  steeds  were 
too    fast. 


amplitudes    cf    all  the    wavenumber    3   features 


63 


V.       CCNCIUSICNS     fiND    SUGGESTIONS    FOP.    FUTURE 
BESrHHrr-H 


lE   this    s"udy   cn€    three-day   and  two   fiva-day    NOGAPS    SOOmb 
forecasts    w€"€    spectrally    verified.  The    cases    were    chosen 

to  reflect  -hree  differer.-  plane-ary  a^a  long  wave  situ- 
a-icns.  In  Case  1  the  5CCmb  patiein  is  characterized  by  a 
vigcrcus  thr':-e  wave  pattern.  In  Case  2  a  forana-.icn  cf  a 
blocking  big!!  and  a  shift  tc  a  low  index  flew  dominate  the 
503[nt  pattarn  while  Case  3  describes  the  breakdown  of  a 
block   and    a    shift  to    a    high   index    flow. 

Two  a.ialvEis  techniques  were  used  to  study  the  observed 
and  forecast  spectral  statistics.  The  amplitudes  of  each 
wavenumber  i";  the  analysis  and  the  forscast  are  compared  at 
key  periods  during  the  forecast  and  Hovmoller  diagrams  are 
prepared  for  key  wav=numker  groupings  following  Baumhefnsr 
and  Ecwney  (1S78).  ^aver.umber  components  were  grouped  into 
planetary  (waves  1-3)  long  (waves  5-7)  and  medium  waves 
(waves   8-12)     and    displayed   on   trough-ridge    diagrams. 

In   all    three      cases    the    most    serious      errors    were    concen- 
trated     m    the      planetary      waves.  The    model      consistently 
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failsd  t.0  corr9c-ly  predict  large  aniplituda  chang'es  in  -hs 
individual  planetary  wavenumbar  components.  NOGAPS  bc-h 
predictai  largs  amplitude  changes  which  did  not  occur  in  the 
observat-ions  and  failed  to  predict  large  amplitude  changes 
which  did  cccur  in  the  analysis.  In  Case  1  the  model  fore- 
cast an  erroneously  large  amplitude  increase  in  wavenumbei  1 
while  in  CaS=  2  NOGAPS  smoothed  the  wavenumber  2  ampli-ude 
which  had  increased  dramatically  in  the  analysis.  Based  on 
these  three  cases  a  systematic  error  can  not  be  established. 
The  erroneous  amplitude  shifts  caused  differences  in  the 
waveruniber  signatures  and  made  estimates  of  phase  speed 
errors  difficult.  Where  dir  ^'ct  comparisons  could  be  made, 
planetary    v»av€    phase    speeds  were    accurately   forecast. 

Th~  Icng  wave  group  (waves  5-7)  was  the  most:  accurately 
forecasc  by  NCGAES.  Cverall  amplitude  and  phase  speed  char- 
acteristics were  well  predicted  although  amplitudes  were 
slightly  siccthed  by  the  ncdel  and  phase  speeds  tended  to  be 
slightly  tec  fast  in  the  regions  of  increased  long  wave 
activity    such    as    in    Case    2.  In    other  regions    of    both    slow 

and  fast  movement  the  model  showed  a  high  degree  of  phase 
speed    accuracy. 
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Ths  medium  wav<=s  (6-12)  verificazion  exhibited  soii^  sys- 
tematic srrors.  NOGAPS  rended  to  smooth  the  amplitudes  cf 
the  arediai  waves  to  the  pcint  vrhera  key  features  were  often 
lost.  Also  phase  speeds  cf  the  medium  waves  were  generally 
too  fast  ifi  the  model.  Specifically,  examples  of  fast  phase 
speeds  in  Case  1  aire  in  the  region  from  OE  to  120E,  in  Case 
2    frcni    130K    tc    30W   and    in    Case   3    from    20E    to    160E. 

Several  sucgestions  for  future  research  follow  from  this 
sxudy.  First,  the  number  of  forecasts  spectrally  verified 
should  be  increased.  New  that  the  spectral  verification 
software  has  teen  prepared,  the  verification  could  be  dene 
routinely.  The  coal  of  the  increased  number  of  cases  would 
be    tc   document    systematic    spectral    errors    in    NOGAPS. 

A  second  area  that  should  be  investigated  is  -^he  rela- 
tionship t€twe€n  -:he  strength  cf  the  circumpolar  vortex  and 
the      forecast      phase      speed        errors.  Variations      in      the 

wavenumber  zero  component  can  be  found  by  comparing  gradi- 
ents cf  the  m^an  height  of  the  forecast  to  the  ODservation. 
The  s-^rength  cf  the  circumpolar  vortsx  may  be  a  factor  in 
the    phase    speeds    of    the    shorter    synoptic    waves. 


66 


A  third  suggestion  involves  tha  quss-ion  cf  "h^  cptiiEui 
north-scuth      ir.-t€rval   for      verification.  This    study      us-d 

eight  degress  of  latitude  for  Case  1  and  five  d=greas  for 
Cases  2  arc  3.  Further  studies  should  examine  -ha  effec-s 
cf  narrc-wer  or  wider  latitude  bands  in  the  analysis.  Simi- 
larly the  effects  of  latitude  variations  on  wavenumbers 
could  be  detcrsined  by  using  data  from  different,  lati-adas 
(60N,   SON,    e-c)  . 

This  study  illustrated  the  utility  of  diagncs-ic  verifi- 
ca-icn  t.echr.iqu<?s  for  study  of  operational  numerical  pre- 
diction medals.  Future  improvemen*:  in  N0GAP3  forecast 
skill,  p  arricularl]'  in  the  medium  range,  is  dependen-  in 
part  en  the  ccntinuir.g  study  of  the  modal  to  isolate  its 
strength    and    weakness  as. 
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